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We apply the nonperturbative functional renormalization group (NP-FRG) in the superfield for- 
malism that we have developed in the preceding paper to study long-standing issues concerning the 
critical behavior of the random field Ising model. Through the introduction of an appropriate reg- 
ulator and a supersymmetry-compatible nonperturbative approximation, we are able to follow the 
supersymmetry, more specifically the superrotational invariance first unveiled by Parisi and Sourlas 
[Phys. Rev. Lett. 43, 744 (1979)], and its spontaneous breaking along the RG flow. Breaking occurs 
below a critical dimension <1dr — 5.1, and the supersymmetry-broken fixed point that controls the 
critical behavior then leads to a breakdown of the "dimensional reduction" property. We solve the 
NP-FRG flow equations numerically and determine the critical exponents as a function of dimension 
down to d < 3, with a good agreement in d — 3 and d — 4 with the existing numerical estimates. 
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I. INTRODUCTION 

Quenched disorder induces an extrinsic inhomogeneity 
in otherwise translationally invariant pure systems. As 
a result, the equilibrium and out-of-equilibrium physics 
of disordered systems may be influenced by rare collec- 
tive events, such as "avalanches", statistically unlikely 
regions, such as "droplets" or "Griffiths regions", and 
by the proliferation of "metastable states". Our objec- 
tive is to develop a theory describing the long-distance 
physics of such systems through a nonperturbative func- 
tional renormalization group (NP-FRG) method and our 
first focus is the equilibrium behavior of the random-field 
model. We showed in previous papers, which we refer to 
as li and 112, that the effect of avalanches and droplets 
can be captured by an approach based on the RG flow 
of the cumulants of the renormalized disorder, provided 
that the full functional dependence of the latter is ac- 
counted fori£ A functional RG is therefore required to 
let the singular behavior due to rare events or regions 
emerge along the flow, as also shown in the case of the 
equilibrium and forced behavior of manifolds in a random 
environment 

The issue of metastable states, which, at zero tem- 
perature where the concept is well defined, refers to the 
presence of many minima of the microscopic hamiltonian 
(bare action) not simply related by symmetry transfor- 
mations, is a recurring conundrum in theories of disor- 
dered systems. In the case of the random-field Ising 
model (RFIM) under study, it has a clear manifesta- 
tion. The critical behavior of the model being controlled 
by a zero-temperature fixed point the long-distance 
physics can be described through the properties of the 
ground state which, due to the presence of the random 
field, is obtained as the solution of a stochastic field 
equation. By standard field-theoretic manipulations, this 
leads to a theory expressed in terms of superfields. Parisi 
and Sourlas 15 showed that a supersymmetry of the the- 



ory, more specifically the invariance under rotations of 
the underlying superspace, implies the property of di- 
mensional reduction, according to which the critical be- 
havior of the RFIM in d dimensions is identical to that 
of the pure Ising model in d — 2 dimensions. The prop- 
erty however has been proven to be wrong in low enough 
dimension ! 16 ' 17 At the same time, it has been understood 
that the superfield construction breaks down because of 
the presence of many solutions of the stochastic field 
equation^ However, in the Parisi-Sourlas formalism, it 
is not possible to disentangle breaking of superrotational 
invariance and collapse of the formalism due to the ap- 
pearance of multiple solutions. 

In the companion paper, henceforth referred to as pa- 
per III , 1 ^ we have shown how to resolve the above conun- 
drum and to combine an extended superfield approach 
with the NP-FRG formalism. In particular, ground-state 
selection can be achieved by adding a weighting factor 
involving an auxiliary temperature and letting at the 
end of the manipulations the auxiliary temperature go to 
zero in the exact NP-FRG equations for the cumulants 
of the renormalized disorder. The resulting property of 
"Grassmannian ultralocality" then allows one to specifi- 
cally investigate supersymmetry (superrotational invari- 
ance) and its spontaneous breaking along the RG flow. 
Such an investigation is the purpose of the present paper. 

The outline of the article is as follows. In Sec. HH we 
briefly summarize the main steps of the NP-FRG in a su- 
perfield formalism that we have developed in paper III. 
This allows us to recall definitions, notations, and RG 
equations that will be used in the present article. In 
particular, we stress two important, and distinct, formal 
properties of the superfield theory: "Grassmannian ul- 
tralocality" and "superrotational invariance". We also 
consider the RG flow of the Ward-Takahashi identities 
associated with the latter and the consequences for the 
choice of the infrared regulator. In Sec. IIII1 we show 
through our NP-FRG formalism that the superrotational 
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invariance nonperturbatively leads to dimensional reduc- 
tion. We conclude the section by building a scenario 
for a spontaneous breaking of the superrotational invari- 
ance along the RG flow that is based on our previous 
results on the appearance of a "cusp" in the functional 
dependence of the cumulants of the renormalized ran- 
dom field^ ' 2 ' 20 i 21 and we propose a continuation of the 
NP-FRG flow equations when spontaneous breaking has 
taken place. 

Sec. |IV] is devoted to the development of a 
supersymmetry-compatible nonperturbative approxima- 
tion scheme for the exact NP-FRG equations. It relies 
on combined truncations of the cumulant expansion and 
the expansion in spatial derivatives of the field. We also 
provide details on the numerical resolution. In Sec.lVlwc 
present the results. We show in particular that break- 
down of the dimensional reduction predictions for the 
critical exponents of the RFIM occurs below a critical 
dimension don — 5.1. We compute the critical expo- 
nents as a function of dimension down to d — 3 and we 
find good agreement with the best available estimates in 
d = 3 and d = 4. Finally, we show that scaling is de- 
scribed by three independent exponents, contrary to a 
proposed conjecture^ - — 

A short account of this work has appeared in Ref . [25[ . 



II. NP-FRG IN THE SUPERFIELD 
FORMALISM FOR THE RFIM 

A. Summary 

We start by briefly recalling the main features of the 
formalism presented in the preceding article, with the 
associated definitions and notations. 

In paper III^ we have developed a NP-FRG theory 
for describing the equilibrium long-distance physics of 
the RFIM that is based on an extension of the Parisi- 
Sourlas 15 supersymmetric formalism. The latter relies on 
the fact that the critical behavior of the model is dom- 
inated by disorder-induced fluctuations, thermal fluctu- 
ations being subdominant, and can therefore be studied 
by looking directly at zero temperature. The equilibrium 
properties are then described by the ground state which 
is solution of the following stochastic field equation 

SS[tp; h] 



5<p(x) 



J(x) 



(1) 



where we have added an external source (a magnetic 
field) J conjugate to the tp field and the action S[ip; h] 
is given by 

S = J x {l VM*)? + Ub(v(x)) ~ hfrMx)}, (2) 

where J x = J d d x, U B {v) = (r/2)ip 2 + (u/4\)(p 4 , with 
h(x) a random (magnetic) field sampled from a Gaus- 
sian distribution of zero mean and variance h{x)h(y) = 
A B SW(x-y). 



Two key ingredients of our extension of the Parisi- 
Sourlas formalism are: 

(1) the need to consider multiple copies (or replicas) of 
the original system with the same disorder, each copy be- 
ing coupled to a different applied source, in order to gen- 
erate cumulants of the renormalized disorder with their 
full functional dependence, thereby allowing for the emer- 
gence of a nonanalytic behavior in the field arguments, 

(2) the introduction of a weighting factor involving an 
auxiliary temperature /3 _1 to the solutions of the stochas- 
tic field equation, so that when /3 _1 appraches only the 
ground state contributes to the generating functional. 

By using standard field-theoretical techniques^ one 
then ends up with a superfield theory for multiple copies 
in a curved superspace, with 



zW[{J a }] = exp(wW[{J }]) = exp(X>f Vj) 

a=l 

P n / rip 

= / nS*aexp( -S {P) [{$ a }] +J2 / Mx)<$> a (3 
J a=l V a=l J * 

whcrc the multicopy action is given by 

-(dMx)) 2 + U B (<S> a (x)) 



(3) 



a=\ J * 

-tEE// *aA^0l)*a 2 (x,9 2 
1 ai =la 2 = l J x J §-ll2 



(4) 



In the above equations, we have introduced a superspace 
(coordinates x) comprising the d-dimensional Euclidean 
space (coordinates x) and a 2-dimensional Grassmannian 
space (anticommuting coordinates = {9, 0}); the metric 
is flat in the Euclidean sector and curved in the Grass- 
mannian one with the curvature proportional to (3. (For 
instance, the integral over superspace is defined as J = 

f x J g = J d d x J f(l+p06)d6d6.) We have also definedsu- 

perfields $ a (x) = Lp a (x)+9ip a (x)+'ipa(x)9+89(i>a{x), with 
one auxiliary bosonic ( "response" ) field <j) a and two auxil- 
iary fermionic ("ghost") fields and ip a , and associated 
supersources J a {x) = J a (x) + 9K a (x)+K a (x)9+99J a (x). 
The action in Eq. (j4]) is invariant under a large group of 
(bosonic and fermionic) symmetries. 

When expanded in increasing number of sums over 
copies, the generating functional of the connected 
Green's functions, [{J a }}, gives access to the cu- 

mulants of the random generating functional [J] . 

We have next applied the NP-FRG formalism to this 
superfield theory. This proceeds by first adding an in- 
frared (IR) regulator that enforces a progressive account 
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of the fluctuations to the bare action, 

=^E / / ^a{x l ,e)R k {\x l -x 2 \)^ a {x 2 ,6) 

L a J xix 2 J0_ 

+ ®a 1 (Xl)Rk{\xi-X 2 \)$ a2 (x 2 ), 

where the two cutoff functions R k and R k are related 
through 

R k (q 2 ) = -^d q2 R k (q 2 ), (6) 

with q the Euclidean momentum, A k the strength of the 
renormalized random field, and Z k the field renormaliza- 
tion constant. This ensures that all symmetries of the 
theory are satisfied. This includes the superrotational 
invariance found when the theory is restricted to a single 
copy and to /3 = 0: it then corresponds to A k = A^Z^, 
a property that is valid at the microscopic (UV) scale A 
(see also below). 

One next introduces the effective average actio n 27 i 28 

if 5 [{*«}] = -W ( f ] [{Ja}] + EaLMaDMx) " 

ASjj/ 3 ' [{$ a }], whose dependence on the IR cutoff k is 
governed by an exact renormalization-group equation 
(ERGE). The effective average action can also be ex- 
panded in increasing number of sums over copies, each 
term of the expansion being then related through the 
Legendre transform to the cumulants of [J] . We re- 
fer to the pth order term of the expansion of [{&a}} as 
the pth "cumulant of the renormalized disorder" . These 
expansions in increasing number of sums over copies lead 
to systematic algebraic manipulations that allow one to 
derive from the ERGE for r^[{$ Q }] a hierarchy of cou- 
pled ERGE's for the cumulants of the renormalized dis- 
order. 

We have unveiled an important property of the random 
generating functional [J] : the latter satisfies a spe- 
cific dependence on the Grassmann coordinates, which 
we have called "Grassmannian ultralocality" , if and only 
if a unique solution of the stochastic field equation is 
included in its computation. This translates into an "ul- 
tralocal" property of the cumulants of the renormalized 
disorder. "Grassmannian ultralocality" becomes a prop- 
erty of the superfield theory when /3 — > oo; it is also 
asymptotically found for finite /3 when k — > (after going 
to dimensionless quantities). It then reflects the desired 
ground-state dominance. When this property is satisfied 
{e.g. by setting =0), the ERGE's for the cumulants 
simplify and only involve " ultralocal" quantities that can 
be evaluated for physical fields $> a {x) = <fia(x), i.e. for 
superfields that are uniform in the Grassmann subspace. 

For illustration, we give below the ERGE's for the first 
two cumulants under the property of "Grassmannian ul- 



tralocality" , which will be needed in the following: 
d t T k i [fa] = 

Jx 2 x 3 

(7) 

and 

d t T k2 <t> 2 ] = 

\dtl {- *W< [<h\ r£i.,* 4 [0i, faM + 

4 (8) 

Pk;x 3Xi [01, 0l] r fe2;i 3 x 4 ,. [01, 02] + [01, 02] 

x ( r Si„,x4 [0i>0 2 ] -Rkixsxt) +perm(12)}, 

where perm(12) denotes the expression obtained by per- 
muting 0i and 02, t = log(fc/A), dt is a short-hand no- 
tation to indicate a derivative acting only on the cutoff 
functions (i.e., dt = dtR k 6/SR k + dtR k 5/SR k ), and the 
propagators P k and P k (denoted by Pj^ and P^ in paper 
III 19 ) are defined as 

P k [4>] = (r k 2) M + R k y 1 (9) 

and 

Pk[fa^2] = PM(r k 1 > l\fa,h]-R k )Pk[h], (io) 

and superscripts indicate functional differentiation with 
respect to the field arguments. 

Generically, the flow of T kp [<p±, ...,(f> p ] involves three 
types of quantities: the propagators P k and P k , second 
functional derivatives of T kp in which all the arguments 
are different, and second functional derivatives of T fc ( p+1 ) 
with two of their arguments equal to each other. A graph- 
ical representation of the hierarchy of ERGE's is provided 
is Appendix C of the companion paper III. 

In the present paper we will focus on these flow equa- 
tions, which correspond to the (3 — ¥ oo limit and de- 
scribe the equilibrium properties associated with the 
ground state (recall that the "bath" temperature is also 
at T = 0). We will briefly comment on the RG flow 
of corrections to "Grassmannian ultralocality" when /3 is 
large but finite. 

Finally, it is worth stressing that to obtain the flow 
equation for T kp [cj)i, <p p ] with its full functional depen- 
dence on the p field arguments, one needs to consider at 
least p copies. If one considers less than p copies, the 
equation necessarily involves T kp and its derivatives in 
which several of the arguments are equal. As we will 
show later, breaking of supersymmetry and breakdown 
of dimensional reduction are precisely related to the an- 
alyticity properties of the fimctionals when taking the 
limit of equal arguments. Formally, the whole hierarchy 
of flow equations for the cumulants can thus be obtained 
by considering an arbitrary large number of copies. 
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B. Superrotational invariance, Ward-Takahashi 
identities and RG flow 



When "Grassmannian ultralocality" is satisfied, the 
curvature /3 disappears from the flow equations [see 
Eqs. (|7|8p ]. Formally, the latter are then the same as 
those obtained from considering a flat Grassmann sub- 
space. We have shown in paper III that in this case, 
and provided one restricts the supersources such that 
one effectively recovers a one-copy theory, the theory is 
invariant under "superrotations" that mix the Euclidean 
and Grassmannian sectors of the superspace. The as- 
sociated generators are Q M = —x^de + (2/Ab)05 m and 
Qfi = x^dg + (2/ Ab)08^. As any linearly realized contin- 
uous symmetry, superrotational invariance leads to a set 
of Ward-Takahashi (WT) identities for the one-particle 
irreducible (1PI) generating functional, i.e. at scale k 
the effective average action. For a flat superspace and a 
restriction to one copy, the WT identity for the effective 
average action reads 

[Hx)Q^im=o (ii) 

J x_ 

and similarly with Q^. 

One can now check that the above WT identities are 
a priori stable under the RG flow, i.e. 



J x 



0, 



(12) 



where Q x generically indicates a component Q M and, 
since there is no curvature, J = j d d x J dOdO. To prove 

Eq. (|12|) . we rewrite the ERGE for the effective average 
action when the superfield theory is restricted to one copy 
and is considered in a flat superspace (j3 = 0): 



where the regulator function TZk;x,,x = ^§ 1 ,e 2 Rk(\xi — 

x 2 \)+Rk{\xi-x 2 \) with S Si S2 = (0i-0 2 )(<?i-0 2 ) and the 
modified (full) propagator Vu is the inverse in the sense 

of operators of V\ ' + T^k- (Inserting the "multifocal" 
expansion in the Grassmann coordinates and using the 
assumption of "Grassmannian ultralocality" then leads 
to ERGE's for the cumulants which are just Eqs. (I7I8[) 
and their extension to higher orders with all copy super- 
fields taken as equal in both sides of the equations.) 

The flow of the left-hand side of Eq. (|12l) can then be 
expressed as 



i 



(14) 



where we have used that Q £i does not act on Vk-x^- 
Assuming that the WT identity, Eq. (ITT1) . is satisfied 
down to the scale k (we study its further evolution), one 
has 

(15) 

which is obtained by differentiating twice Eq. (fTTj) . Af- 
ter an integration by parts and a relabel of the dummy 
variables x 2 , x 3 , one then finds 

= -\*l T i%* 2 [ $ ] (q £2 + 2*3) w, Vl m . 

Jx 2 x 3 

Provided that one chooses the infrared cutoff function 
such that 



^Q Xl + Qx 2 ) T^k-x ± x 2 — 0, 
-(2) frRl u„ -r( 2 ) 



(17) 



one can replace r^ £ J$] by r^jS] + TZ k ,^ 2 in 
Eq. p^|) . By using the fact that the latter is equal to 
the inverse modified propagator (P^ 1 ^})^^, one can 
rewrite the right-hand side of Eq. (fT5|) (up to the trivial 
factor of —1/2) as 



•t 1 y , 



d t T k {{<£ a }} = \ [ {Wk^a,) (13) = q f 

J Xr, 



J x n x* J x„ 



(18) 



This last expression is easily shown to be identically zero, 
so that one finally obtains that Eq. (|12l) is satisfied at 
scale k. The same reasoning can be repeated with the 
identity associated with Q x . 

Therefore, if one starts with an initial condition that 
is superrotationally invariant, which can be enforced by 
suppressing fluctuations (see section II- A and paper III), 
the above derivation proves that the WT identities as- 
sociated with the superrotational invariance are a pri- 
ori preserved along the RG flow if the regulator satisfies 
Eq. (fTTj) . The latter condition can be reexpressed as 



2 ~ 
+ -t— d^Rk^X-y - x 2 \ 



(Xl,t - X2fj.)Rk(\xi - x 2 \) 



(19) 



where we have used the fact that R k and R k are transla- 
tionally invariant functions of the Euclidean coordinates, 
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with therefore di^Rk — —d-^Rk and similarly for Rk- 
Going to Fourier space, this implies that Rk and Rk are 
related through 



Rk{q 2 



-A B d q 2R k (q 2 ), 



(20) 



which, as already stated, corresponds to Eq. © with 
Afe = Afl-Zfc, a condition that we choose to enforce at 
the beginning of the flow when k — A. Such a choice of 
cutoff then avoids an explicit breaking of superrotational 



mvariance. 



III. SUPERROTATIONAL INVARIANCE AND 
ITS SPONTANEOUS BREAKING 

A. SUSY and dimensional reduction 

An important property of the theory is that the super- 
rotational invariance (which for simplicity will often be 
denoted by the acronym SUSY, for supersymmetry, in 
the following) leads to dimensional reduction. We con- 
sider here the case where "Grassmannian ultralocality" 
is satisfied so that the curvature /3 drops out of the flow 
equations as in Eqs. (|7I8I) . and we assume that SUSY is 
indeed obeyed when the theory is restricted to one copy 
(see above) . We showed in paper that SUSY implies 
nontrivial WT identities relating cumulants of different 
orders. More specifically, the following WT identity at 
the scale k will be needed below: 



A 



M\p( 2 ) 



(21) 



To prove the dimensional reduction property, it is then 
useful to single out a 2-dimensional subspace of the d- 
dimensional Euclidean space. We define x = (y, z) with 
y € R d ~ 2 and z € R 2 and we consider superfields that 
are uniform in the 2-dimensional Grassmannian and the 
2-dimensional Euclidean subspaces, i.e. §{x) — <fi(y). 
For such fields, the WT identity in Eq. (f2"Tj) implies that 



.(ii) 

fe2;y 1,1/2 



[P 2 



= A B d„2T 



(2) 

kl;yiy 2 



(22) 



where p is the 2-dimensional momentum obtained from 
a spatial Fourier transform over z [by using transla- 

(2) 

tional invariance we have as usual defined F,J {pi,P2j — 
(27r) 2 8^( Pl +p 2 )T [ ^{p 2 ) and similarly for T^\ u p 2 )]. 
The same type of relation holds between —Rk^t^ip 2 ) 
and Rk lVl ,y 2 (p 2 ). 

With the help of the above relations, the ERGE for 



the first cumulant r^i [<j>], Eq. 0, can be reexpressed as 

id-2„. [ d2 P 



d t r kl [<f>] = -(^-)d t 



i/i 



*Pk Wl yAp 2 ,4>\d p 2 ^ { kly iy2 [p 2 



(2tt) 2 

^Rk( P 2 : 



x a, 



47T 

(2)r 2 M , 3 fJ2' 



log r&V,d + W 



(23) 

(2) ^ 

where, in the last expression, T k ^ and Rk are functions 
of p 2 but operators in the (d — 2)-dimensional Euclidean 
space spanned by the coordinate y (since <fi(y) is not uni- 
form) . The integral over p 2 is easily performed. Choosing 
an infrared cutoff function that becomes independent of k 
when its argument p 2 is at the UV scale, it only remains: 



d d ~ 2 y 



x [log {r^[p 2 = 0, 0] + R k (p 2 = 0)) ] yy . 

(24) 

Up to the trivial factor As/(47r) and with the identifica- 
tions 



T kl [cj>} = T k m,T {2 }[p 2 = 0,0] = T k A >[<f>],R k (p* =Q)=R k , 

(25) 

Eq. (f24|) coincides with the ERGE for the effective av- 
erage action of a standard scalar cf) 4 theory 28 in di- 
mension d — 2. This provides another nonperturbative 
demonstration 2 ^— of the property of dimensional reduc- 
tion for supersymmetric scalar field theories, first put for- 
ward by Parisi and Sourlas. 15 



(2)r 



B. Spontaneous SUSY breaking 

We have shown within a nonperturbative implemen- 
tation of the RG that the supersymmetry of the theory, 
and more specifically the superrotational invariance for 
one copy (SUSY), imply the property of dimensional re- 
duction. As one knows that dimensional reduction does 
not hold in low enough dimension, what then goes wrong 
in the formalism? 

The answer is that some (super)symmetries or iden- 
tities must be spontaneously broken along the RG flow. 
We make the assumption, which will be supported by ac- 
tual calculations, that only the superrotational invariance 
(for the theory restricted to a single copy) may be broken 
along the RG flow, all other (super)symmetries and prop- 
erties, most significantly the "Grassmannian ultralocal- 
ity" encoding ground-state dominance (when /3 —¥ oo), 
remaining unaltered. From the above proof, it follows 
that failure of dimensional reduction implies that the WT 
identity between and breaks down: a singular- 
ity must occur along the flow, which invalidates the WT 
identity; the equality in Eq. (|12[) and some step in its 
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derivation must go wrong at some scale k. From our pre- 
vious workf^ i 20 ' 21 we anticipate that this arises due to 
the appearance of a strong enough nonanalytic depen- 
dence of the second cumulant of the renormalized ran- 
dom held r^^i, 4>2\ in its field arguments when the 
latter two, and 4>2, become equal. (One should keep 
in mind that the cumulants are invariant under permu- 
tations of their arguments; consequently, rj^ is an even 
function(al) in — 2 -) As a result, some higher-order 
derivative (i.e., a 1PI vertex deriving from rl^ ) blows 
up and the whole hierarchy of coupled ERGE's and WT 
identities for one copy ceases to be valid. 

In order to envisage the possible scenario for sponta- 
neous SUSY breaking and breakdown of dimensional re- 
duction, it is instructive to study the structure of the ex- 
act hierarchy of coupled RG flow equations, whose lowest- 
order examples are given in Eqs. (|7l8p . The following ex- 
position is not meant to be rigorous, but only to provide 
some heuristic arguments. 

Consider first the hierarchy of ERGE's for the proper 
vertices evaluated for uniform fields (in Euclidean space) 
that is obtained by repeated functional differentiation of 
the ERGE for the effective average action. For this part 
of the reasoning, the case of a simple scalar field theory 
is sufficiently illustrative. The ERGE for the effective 
average action then reads 

W0] = I [ d t R k (q 2 ) (ri 2) [0] + RX 1 . (26) 
* Jq v ' g,-g 

After differentiating this equation s times, it is easily re- 
alized that for s > 2 the ERGE for any generic 1PI vertex 
rl depends on all lower-order proper vertices (of order 

> 2) as well as on some F^ s+1 - ) 's and F^ s+2 - ) 's. Two cases 
then arise: when s > 5, the flow equation is linear in 

(s) 

T k itself, whereas in the other case when 4 > s > 2, the 

is) 

flow equation may be nonlinear in T k itself. The same 
structure applies to the present theory. A consequence is 

(s) 

that, if no other 1PI vertices diverge first, T) can only 
diverge at the end of the flow (k — > 0) when s > 5 (due 
to the linearity of the corresponding ERGE) whereas a 
divergence may occur at a finite scale when s < 4. 

To get more insight, we have to look more thoroughly 
into the structure of the RG equations and consider the 
1PI vertices associated with the cumulants separately [see 
Eqs. (|7I8I) for illustration]. We further make the following 
plausible assumptions (in order to have a renormalizable 
theory) : 

(i) The 1-copy 1PI vertices rl* always stay finite, as do 

the cumulants of the renormalized random field ^^p" 
and their first-order derivatives. The latter may however 
be discontinuous; the condition on the first derivatives 
amounts to assuming that no "supercusp"— stronger 
than linear appears in the cumulants when two argu- 
ments become equal. 

(ii) Since we expect a nonanalytic behavior occur- 
ing as the copy fields become equal in the reduction 



to an effective 1-copy system, we also anticipate that 
derivatives with respect to symmetric combinations of 
the field arguments are always bounded (e.g., \(d(j>i + 

d</ , 2)rj.p lS2 '" Sp) [</!>i! 4>2, ■•■<MUI < +oo, where the deriva- 
tive is evaluated for all replica fields equal to tp). 

A "cusp"— in r? 1 '" 1 -' with p > 3 implies that one of 
its second derivatives, which means a proper vertex of 
order at least s = 5, blows up. Similarly, for having a 
"subcusp"— in rj^ , at least one of its third derivatives, 
hence a proper vertex of order at least s — 5, should 
blow up. As we have seen above that the associated flow 
equations are linear, we conclude that the divergence of 
such proper vertices (with s > 5) cannot be the trig- 
gering event for a singularity at a finite scale. On the 
other hand, due to the nonlinearity of the ERGE's for 
the second derivatives of T^P, a "cusp" may appear in 
the latter at a finite scale. 

To summarize this discussion: one expects that break- 
down of dimensional reduction requires the presence of a 
cusp in the field dependence of rj^, cusp that should ap- 
pear at a finite scale, which by analogy with what occurs 
for a manifold in a random environmen t 10 ' 33 we call the 
"Larkin" scale, during the RG flow. By contrast, weaker 
nonanalyticities in rj^ as well as nonanalyticities in 
higher-order cumulants (if no cusp has appeared in rl^) 
can only appear at the fixed point, in the limit k — > 0, 
thereby preserving the dimensional-reduction property. 
Indeed, the derivation leading to Eq. (l24l) then remains 
valid. 



C. Continuing the NP-FRG flow with broken 
SUSY 

Another nontrivial question is then raised: if superro- 
tational invariance (SUSY) is spontaneously broken, how 
can one continue the RG flow for the effective average 
action? The answer, which again should be verified in 
calculations, is as follows: 

(i) Keeping in mind that superrotational invariance is 
explicitly broken in the presence of multiple copies (see 
section [TTJ) and that it is only recovered in the process 
of reducing the whole problem to a single-copy system, 
spontaneous breaking may then be described by reject- 
ing any implicit assumption of analyticity of the field 
dependences in the latter process. One should therefore 
restrict the hierarchy of ERGE's to those equations for 
cumulants that are considered for generic (nonequal) field 
arguments, so that a putative nonanalytic dependence in 
these arguments can freely emerge along the RG flow.— 

(ii) One assumes that except for the superrotational 
invariance, all of the properties and symmetries of the 
effective average action remain valid. In particular, 
this applies to the property of "Grassmannian ultralo- 
cality" of the random generating functional and to its 
consequences for the cumulants of the renormalized dis- 
order. As discussed in the preceding paper, this property 
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is obeyed when taking the limit = in the ERGE's 
for the cumulants, which leads to Eqs. (|7I8[) and their ex- 
tensions to higher orders. It keeps track of the fact that 
only the ground state is considered for each copy. The 
validity of the property of "Grassmannian ultralocality" 
is therefore distinct from that of superrotational invari- 
ance, a distinction that cannot be made in the original 
Parisi-Sourlas formalism. 

(iii) The cutoff functions satisfy the relation in Eq. ([6]) 
which contains as a special case the explicit superrota- 
tional invariance when A k = A^Z/j. In practice, the 
typical strength of the renormalized random field Aft is 
defined from Vji2 (q 2 = 0) evaluated for some specific 
uniform field configuration and the field renormalization 
constant Z k is obtained from T$(q 2 = 0) also evalu- 
ated for some uniform field configuration. Finally, the 

cusp in rj^j. X2 {4>i,4>2), which is associated with spon- 
taneous SUSY breaking, must be stable upon further 
evolution with the infrared scale k and must not gen- 
erate "supercusps" . From the structure of the ERGE 
for Y^) xi (cfrijfc), h is expected that a "linear cusp"— 
satisfies these requirements and provides a mechanism 
for dimensional reduction if it persists at the fixed point 
(see Appendix |A"t . 



to describe the physics of the RFIM. The first order of 
the derivative expansion leads to T^(q 2 ;<j)) of the form 
U'k{4>) + Zk((f>)q 2 , with Uk(<t>) the "effective average po- 
tential" ( "Gibbs free energy" for a magnetic system) and 
Zk(4>) a field renormalization function. When the WT 
identity is satisfied, this in turn requires that 



T^\q 2 -^(j>)=^ B Z k {4>), 



(27) 



which corresponds to a local approximation for the sec- 
ond cumulant. It is easy to generalize this reasoning by 
taking into account the whole set of WT identities for the 
proper vertices obtained from the cumulants: a SUSY- 
compatible approximation at order n consists of taking 
Ffci at the order n of the derivative expansion, at 
the order n— 1 of the derivative expansion, T).r n+ i-\ in 
the local approximation, and all higher-order cumulants 
equal to zeroi2£i This provides a scheme of successive ap- 
proximations that in principle can be used to check the 
robustness of the results and if necessary improve them. 

The minimal truncation that can already describe the 
long-distance physics of the RFIM and does not explicitly 
break SUSY is then the following: 
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u k (4>(x)) + -z k ^(x))(d^(x)Y 



(28) 



IV. APPROXIMATION SCHEME 
A. SUSY-compatible approximation scheme 

Up to this point, many of the previous considera- 
tions remain plausible conjectures. We now provide 
an implementation of the NP-FRG that allows us to 
check their validity. The hierachy of ERGE's of course 
cannot be solved exactly and the next step is to pro- 
vide a SUSY-compatible nonperturbative approximation 
scheme. From our previous work ; 1 ' 2 ' 20 ^ 1 we know that an 
efficient scheme relies on a joint truncation in the deriva- 
tive expansion, which approximates the long-distance be- 
havior of the 1PI vertices,— and in the expansion in cu- 
mulants of the renormalized disorder. The truncations 
however must be combined in a way that does not ex- 
plicitly break the supersymmetry, more precisely the su- 
perrotational invariance (SUSY). To implement this re- 
quirement, we use the WT identities. 

As is clear from the discussion in the preceding sec- 
tion, a minimal truncation must at least include the sec- 
ond renormalized cumulant. When SUSY is not bro- 
ken, the WT identity [see Eq. (|22[) ] imposes that for a 

uniform field (q 2 ;4>, 4>) is given by the derivative of 

T$ (q 2 ; 4>) with respect to q 2 . (A similar relation is sat- 
isfied by the infrared cutoff functions.) Consider now 
the derivative expansion. In the lowest order (LPA for 

(2) 

"local potential approximation"), the derivative of ri, J 

is a field-independent constant, which would force 

to be also field independent: this is clearly too crude 



- k2\ 



(29) 



with the higher-order cumulants set to zero.— 

Inserted in the ERGE's for the cumulants, Eqs. 
and ([5]). the above ansatz provides 3 coupled flow equa- 
tions for the 1-copy effective average potential U k (4>) (or 
its derivative) that describes the thermodynamics of the 
system, the field renormalization function Z k (<p), and 
the 2-copy effective average potential V k (<pi,<p2) from 
which one obtains the second cumulant of the renor- 
malized random field at zero momentum, r^ 1 ^? 2 = 



Q-4>i,h) = ^i^2) = V h 



(11)/ 



b 2 ). SUSY is obeyed 



when A k (4>. 
UV scale k 



A. 



Z k (4>), which is easily satisfied at the 



B. NP-FRG equations in a scaled form 

In order to search for the fixed point that controls the 
critical behavior, the flow equations must be recast in a 
scaled form. This can be done by_ introducing appropriate 
scaling dimensions (see refs. [1,20] and paper III 19 ). Near 
a zero-temperature fixed point ) 12 i 13 one has the following 
scaling dimensions: 



Z k ~ k-\ A k ~ k-^-V, 4> a ~ k^ d - 4+f », 



(30) 



and a renormalized temperature is introduced as T k ~ 
~ k e //3, with 6 and fj related through 9 = 2 + r) — ry. 



More precisely, we define running anomalous dimensions 
T) k and fj k as 



d t \ogZ k = -rjk, 



d t log A k 



-(2/7fe - Vk), 



(31) 



(32) 



where Z k and Afc have been introduced in Sec. IV-D. One 
also has 



J-k. -*- u 



(33) 



The dimensionless counterparts of U k ,V k , A k ,(f> will be 
denoted by lower-case letters, u k ,v k ,S k , (p. 

The resulting equations are 



d t u' k (ip) = -i(d - 2?7 fc + f] k )u' k (ip) + i(d - 4 + ryfeV 



x „£'(¥>)] [z fc (^ - S k , Q (<p)] + l[ d \<p) [ - ~4(^) + 4 >0 (^)] 



+ \(vk m) [4(v)jf d+2) (v) + }, 



(34) 



d- 2 



d t v k (ipi,ip 2 ) = -(d - 4 + 2% - 2i] k )v k (ip ll tp 2 ) + 
-(d- 4 + f}k)(<Ptd Vl + (p2d V2 )v k (if 1 ,(p 2 ) - 2wd| - x 
'lA^i.Vs)^^!,^) - z fe (<ysi)] 2 + 4 d) (^i)wi 20) (^i^2) 



x [5fc,o(vi) _ z fe(^i)] + 



d-2 



4 d " 2 V>fV^)+ 



-m^j 2) ((^i,^ 2 ) - 2 [h{vi,V?) ~ zfc(^i)] 
+ 4 20) (</>i ) ^2)i2 1,d) (<Pi)] + perm(12) | , 

(36) 

where v^ 1 = 2 d+1 TT d ^ 2 T(d/2), a prime denotes a deriva- 
tive with respect to the field (when only one argu- 
ment is present), 8 k (ip 1 ,Lp 2 ) = v[ n) (^i,<^ 2 ), S k ,o{<Pi) = 
6k((fi, (fi), and perm(12) denotes the terms obtained by 
permuting ipi and ip 2 . 

The functions lq(<f), Iquq^i,^), TriqQp), 
™>qt]q 2 {<Pi>¥>2), n q i] q2 (ipi,cp 2 ), h q d \ip), jqt iv) and 
jq^lql (ipi, pi) appearing in the above flow equations are 
"dimensionless threshold functions" which are defined 
from the infrared cutoff functions that satisfy 35 



and 



R k (q 2 ) = Z k k 2 s{q 2 /k 2 ) 



R k (q 2 ) = -A k s'(q 2 /k 2 ), 



(37) 



(38) 



d t z k (tp) = i] k z k (ip) + -^{d - i + rj k )(p z' k (tp) - 2« d | 

l[ d - 2 \<p)4(<p) (d- 2)4 d - 2 V)4(¥>K'(¥>) + 2 ^x 
lf +2 \ v )z' k { v f[z k { v ) - 6 k , Q (v)} +4l¥\<p)z' k (<p)u' k "(<p)[z k (<p) 



SkAv)] - 4 d) (y) P" r - 4(y) 2 + z'k(<p)( z k(<p) - 4,o (y)) 



-2d+l , 2 



- 24(^)^,0 (¥>)] + m i (vWkiv) 2 + -^A (^X'(^)x 

[(d + 2)4(^) - 2^ fo>)] + imi d+2 V)4M[(d + 4)4(^)- 
4^»] - -Amf +A \ V )z k {yf + 2mf +2 \ V )z k {v)u'£ { V )+ 



^Wt^l^W-^oM] - o(% - %)[4'M x 



- (4 + ^)4(^) 2 if ' d+2 V) -44(^)<(^)x 

i M V) + ^4(^) 2 4 d+2 V) + ^4(^)< , (^4 d) (^+ 



with the prime again denoting a derivative. The de- 
pendence of the dimensionless threshold functions on 
the dimensionless fields comes through u'^(ip) and z k ((p). 
Their definitions are given in Appendix [B] The proper- 
ties of these threshold functions have been extensively 
discussed i 28 ' 39 ' 40 They decay rapidly when their argu- 
ments u' k '[ip a ) become large, which, since u k '{ip a ) = 
U k (4> a )/(Z k k 2 ) is the square of a renormalized mass, en- 
sures that only modes with a mass smaller than k con- 
tribute to the flow in Eqs. (|34M36I) . These threshold func- 
tions essentially encode the nonperturbative effects be- 
yond the standard one-loop approximation! 28 ' 39 ' 40 

Note that u k (ip) and z k (tp) are even functions of (p 
and that, due to the Z 2 symmetry and the permuta- 
tional symmetry of the cumulants in their field argu- 
ments, Vk(<pi,<pv) satisfies v k (ip ly (p 2 ) = v k (-(pi,-tp 2 ) = 
Vk(<P2,<Pi) = v k (-(f2,-(pi), and simil arly for jhJjfu <fi ;)- 
Finally, the initial conditions for Eqs. pi!32l34l35l36p at 
the UV scale A (t = 0) are obtained from the bare action 
in Eq. (g]) for ,3 = and $ a = ip. a , i.e., 



z A {ip) = 1 



4! 
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4! 



(^ 2 - «a) 2 , 



(39) 
(40) 



(35) 



the last equality being valid in the (relevant) region where 
the bare potential has a nontrivial minimum; 28 after a 
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trivial rescaling of the fields so that the explicit depen- 
dence on A# disappears, one also has 



(41) 



which implies 5\ (tpi, ^2) — 1- This initial condition triv- 
ially satisfies 5a{<P, <p) — z\(p) which ensures that invari- 
ance under superrotations is obeyed at the beginning of 
the flow. 



C. Properties of the flow equations 

An important property of the present theory occurs in 
the limit (p2 — > <p>\ when 4(</ 3 i, ¥2) is regular enough. For 
convenience, we introduce the variables ip = (ipi + ¥2)/% 
and y — (pi — tp2)/2. The property can be stated as 
follows: If the behavior of 5 k is regular enough when 
y — » 0, i.e., is such that 

s k (ip, y) = 4,0(95) + ^MfoOw 2 + y 2 °(y) ( 42 ) 

with o(y) — > when y — > 0, then the flow of 5 k fi(ip) co- 
incides with that of z k (p). This is precisely the WT 
relation in Eq. (|22p when evaluated at zero momen- 
tum; indeed, in dimensionless form, this equation is just 
4,o(</?) = z k (ip). Within the present ansatz for the effec- 
tive average action, dimensional reduction then follows. 

To prove the equality of the two flow equations, we 
first derive Eq. (|36p with respect to <p>\ and (f2 to obtain 
the RG equation for S k , change variables to <p, y and take 
the limit y — by assuming the above small y behavior 
[Eq. (|42|) ]. The derivation makes use of the properties 
of the dimensionless threhold functions that are summa- 
rized in Appendix [B] This leads to 



rd-2 ' 
- 2vA - 



(rj k - 2r) k )6 kfi ((p) + -(d-4 + i] k )cp6' kfi (tp) 



(d-2). 



2d -3 



u'k (tpWk.oip) + m 4 d+2) ( ( P)<5fe,o( 1 ^) 2 + L k[p] 4,0 - z k ]j, 

(43) 

where L k [tp;S k ,o ~ z k ] is a polynomial of degree two in 
[4,0 — z k (tp)] and its derivatives, with coefficients that 
are function of ip and such that L k [ip; 0] = 0. Subtracting 
from the above flow equation the one for z k in Eq. (|35[) 
leads to an equation for the difference 5 k ,o — z k in the 
form 

dt [4,o (<^) ~ z k (ip)] = M k [p\S k fi - z k ] - (fjk - Vk)zk{p), 

(44) 

where M k [ip;S kt o — z k ] is another polynomial of degree 
two in [4,o (v) — z k ((p)] & n d its derivatives such that 
M k [ip; 0] = 0. [Note that due to their definitions, f) k and 
r) k are equal as soon as 4,o(v) = z k (ip).] Since the initial 



(d-2) 



id), 



condition satisfies 4=a(vi, P2) = 1 and z k =h{ l p) = 1 5 
the only solution to Eq. (|4"4"|) is 4,0 (v) — z k (<p) = at 
all scales, which provides the desired demonstration. By 
inserting this result in Eq. (|3~4")l . one obtains 

dtu' k (<p) = -~(d - Vk)u' k (<p) + -{d - 4 + r) k )(f u k (<p) 



(d -2)v d {l { t 2) (p) u'i' (*>) + l[ d) (ip) z' k (p)}, 



(45) 



which coincides with the flow equation for u' k ((p) in the 
<p A theory without random field in dimension d — 2 (use 
the relation (d — 2)vd = (2it)~ 1 Vd~2 and compare with 
Eq. (4.33) of Ref. [28]). The same exercise can be re- 
peated for the flow of z k ((p), thereby completing the proof 
of dimensional reduction in this nonperturbative approx- 
imation. 

On the other hand, a spontaneous breaking of SUSY 
and of the associated WT identity occurs whenever 
4,2 (<p) diverges and 5 k has a cusp-like singularity in the 
form 



4(</>i, <P2) = 4,o(<p) + 5 k ,i(p)\pi - V?2| 



(46) 



as (p2 —> pi — >• ip. We expect that such a cusp appears 
at a finite "Larkin" scale k^ through the formation of a 
boundary layer as k — > kh and y — > in the form 



d 2 y 5 k (ip,y) 



1 



[k - k L (p)] 2 ™ + y 1 



(47) 



with w an (a "priori unknown) exponent, as can be 
checked by inserting in the flow equations. The bound- 
ary layer provides the appropriate initial conditions to 
further continue the flow below k^ with a cuspy func- 
tion. 

Finally, it is worth comparing the above NP-FRG 
equations with those we have previously derived without 
the superfield formalism 1 '- 2 in a minimal approximation 
based on the same ansatz as in Eqs. (|28I29|) . The two 
main differences are that: 

(i) The present equations are directly considered at 
zero temperature (and zero auxiliary temperature) and 
subdominant terms proportional to T k do not appear in 
the beta functionals. 

(ii) More importantly, since from the previous formal- 
ism we had no insight into the requirements for avoiding 
an explicit breaking of the underlying superrotational in- 
variance at the single-copy level, we had for simplicity set 
the infrared cutoff function R k to zero. As a result, the 
WT identity associated with superrotational invariance 
was violated at all scales. The dimensional-reduction 
property could therefore not be exactly recovered when 
it was expected to be validpii 2 - 



D. Numerical resolution 

We are interested in the long-distance physics of the 
RFIM near its critical point and in the associated (once 
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unstable) fixed point. Such information can be obtained 
from the solution of the above coupled FRG flow equa- 
tions with appropriate boundary conditions. We have 
carried out the three following types of resolution. 

(1) We have investigated the location of the break- 
down of SUSY associated with the appearance of a cusp 
in 5k((fi,ip2)- To this end, we have solved the set of 
coupled flow equations which is obtained when 5k has a 
sufficiently regular (no cusp) behavior as tpi — ¥ 932, 

y —> 0: see Eq. (|4"2"j) . As already shown, the flow of 
<5/c,o(<p) then coincides with that of Zk (93), and the NP- 
FRG equations for u' k ((p) and z k (<p) form a closed set. 
In addition, we have solved the equation for the second 
derivative of 5k with respect to (ip% — 932), <5fc,2(v), which 
depends on u' k ((p) and Zk{f)- A divergence in 5fc,2(y) 
signals the appearance of a cusp in Sk(<fi, <P2)- 

The numerical task of solving this system of three cou- 
pled partial differential equation is reasonable. We dis- 
cretize the problem both in the time and field directions. 
We use finite differences to evaluate the field derivatives 
that appear in the beta functions and use an explicit 
scheme to find the (RG) "time" evolution of the func- 
tions. 

(2) When a cusp is present, solving in their full glory 
the set of coupled partial differential equations for func- 
tions of up to two field arguments is a hard task. We have 
discretized the field variables, using a 2-dimensional grid 
to solve the equation for 5k(ifi, if 2)- For convenience, we 
have switched variables from tpi,tp% to X = 931932/2 and 
y — (ipi — (fi2 )/2 and restricted the grid to a trapezoidal 
region with X > and y > 0. Observe that the flow 
of 5k(<fi,<P2) depends on 8k(<pi,<fi) and 6^2,^2) and 
their derivatives. In this sense, the flow equation for 5k is 
nonlocal since its evolution depends not only on its prop- 
erty in a neighborhood of (931, 932) but also on its values 
at points faraway in the (931,932) plane. It is therefore 
important to choose the trapezoidal grid such that, for 
all point in the grid, the information necessary for com- 
puting the evolution of the function at this point is inside 
the trapezoid. In the larger simulations we had 60 points 
in the X direction and 50 in the y direction. 

We have solved the set of coupled equations with two 
kinds of modus operandi: 

(2-i) First, for selected values of the dimension d we 
have searched the "cuspy" fixed point by dichotomy, 
starting from "cuspy" initial conditions (5k=A,i(%) 7^ 0). 
Depending on the dimension and the regulator param- 
eters, there appears from time to time numerical insta- 
bilities. Such instabilities showed up for instance in the 
(physically not very interesting) large-field region, where 
the threshold functions are very small. To keep the latter 
under control, we have therefore used slightly modified 
flow equations: this is discussed in Appendix [Cl 

(2-ii) Secondly, once a proper "cuspy" fixed point has 
been identified in some given d, we have followed this 
fixed point by continuously decreasing or increasing d. 
This was done by looking directly for the roots of the 
beta functions by means of the Hybrid algorithm, using 



as an initial guess the fixed point found at the previous 
iteration. We have studied the cuspy fixed point down to 
d < 3; it is indeed difficult to go down to lower d because 
the anomalous dimensions r\ and fj become large in low 
dimensions: fj ~ 4 — d (see also Ref. [Hj]). 

(3) When approaching the critical dimension d^R for 
both SUSY and dimensional-reduction breaking from be- 
low, the procedures (2-i) and (2-ii) above become ineffi- 
cient. We have used instead an expansion of the function 
<5fe(93i = x + y, 932 = x — y) in \y\ with, in practice, 

V 2 

6k(x+y,x-y) = ^. (a;) + |y|^ 4 (x) + ^-4,2(a:) + O(|?/| 3 ). 

(48) 

Note that the expansion cannot be applied in the close 
vicinity of dun- Indeed, one then expects the approach 
d — > dp R to proceed nonuniformly in y through a bound- 
ary layer in y 2 /(doR — d). It is easy to check that a 
similar boundary-layer phenomenon is found in the more 
easily handled 1-loop equation for the RFO(iV)M near 
d = 4,— around the critical value Ndr — 18, i.e. when 
N -> 18-. 

To compute the anomalous dimensions, we have de- 
fined the amplitudes Zk and appearing in Eqs. (131I32[) 
by choosing the specific configuration of the fields as 
equal to zero. This amounts to imposing 

d t z k (<p = Q) =0, (49) 

d t 8 k ,a((p) = d t Sk(>p, 95) = 0, (50) 

which fixes r/k and fjk- 

In addition to the anomalous dimensions r\ and fj (ob- 
tained as the fixed-point values 77* and 77*), we have com- 
puted the critical exponent v that controls the escape of 
the RG flow away from the fixed point along the unstable 
(relevant) direction. This was done by diagonalizing the 
flow equations around the fixed point and extracting the 
negative eigenvalue which identifies with — 1/v. 

Finally, we have chosen a dimensionless cutoff function 
s{q 2 /k 2 ) [see Eqs. (|37} and (J3SJ)] of the form 

s(y) = (a + by + cy 2 ) e~ y , (51) 

where y = q 2 /k 2 [not to be confused with (932 — 9?i)/2]. 
This function satisfies the requirements s(y — > 00) — > 
and s(y — > 0) = est > 0. It is a modification of the 
function used in Refs. p7ll28| . The parameters a, b 
and c can be optimized via stability considerations (see 
Refs. 43.44.45]). They are varied to find a region of val- 
ues for which the computed critical exponents are stable 
under changes of the parameters. For more details, see 
Appendix [Cj 

V. RESULTS AND DISCUSSION 
A. Results 

Our main result is a numerical confirmation of the sce- 
nario put forward above: the NP-FRG equations allow- 
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ing us to continuously follow the critical behavior and 
the associated fixed point as a function of dimension d, 
we find that there is a critical dimension dun separating 
a region d > don where SUSY is valid all along the RG 
flow (except possibly right at the fixed point) and dimen- 
sional reduction applies from a region d < doR where 
SUSY is broken at a finite IR scale along the RG flow 
and a breakdown of dimensional reduction takes place. 
By using the procedure (1) detailed in the above subsec- 
tion, we have numerically located this critical dimension 
as djjR ~ 5.1 ± 0.1 (the precise value has a residual de- 
pendence on the chosen cutoff function). Note that the 
value d]jft ~ 5.1 obtained here is consistent with the 
value found in our previous, and somewhat cruder, NP- 
FRG approach of the RFO(Y)M when extrapolating the 
transition line duR(N) down to N = 1 (see Figure 4 of 
paper II^). 

For initial conditions of the RG flow at or near the criti- 
cal point, the second derivative Sk^iv) — 9y5k(f + y, <p — 
y)\ y =o blows up at a finite RG "time" 7j £ — log(fcz,/A) 
for d < don, whereas it stays finite up to the fixed point 
for d > doR^ We illustrate the difference of behavior 
between these two cases in Fig. [TJ for a field configura- 
tion ip — 0. The divergence of the full function Sje^if) 
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FIG. 2: NP-FRG flow of 8 k , 2 (p) = d%8(p,y)\ y= o for d = 5 < 
(Idr- The initial conditions at k = A (i.e., t = 0) for u' k (p) and 
Zk{p) are taken at the fixed-point solution, u' t (p) and z,(p), 
and those for Sk,o(p) and Sk,2(p) are chosen as S k= A,o(p) = 
z*(p) and 8k=A,2(p) — 0. One observes that the divergence 
first takes place for small values of p. 
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FIG. 3: Fixed-point solution 8* t 2(p) for dimensions d > doi?. 
Note the square-root behavior V 'd — don as a function of di- 
mension. 
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1*1 

FIG. 1: NP-FRG flow of 5 fe , 2 (0) in the regime where SUSY is 
valid. The initial conditions at k — A (i.e., t = 0) for u' k (p) 
and z k (p) = 5k,a(p), with p — p 2 /2, are taken at the fixed- 
point solution, u'*(p) and z,(p) [d t u' k (p)\ t — dtz k (p)\* — 0], 
and those for 8 k ,o(p) and 8 k ^(p) are chosen as 8 k= A,o(p) = 
z,(p) and 5fc=A,2(p) = 0. The upper (color online blue) curve 
corresponds to d = 5.2 > dmi and one observes that 8k, 2(0) 
tends to a finite fixed-point value. The lower (color online 
red) curve corresponds to d = 5 < doR shows a divergence at 
a finite RG "time" tL- 

when d < d^R is shown in Fig. [2] (Due to the Z2 sym- 
metry, it is more convenient to represent the functions 
in terms of p = (p 2 /2.) We find, as can be anticipated 
from an analysis of the NP-FRG equations for u' k (ip), 
Skfli^) and Sj s: ,2{ l p)i that the latter stays finite at d = doR. 
and that its fixed-point value for d d~jj R behaves as 

S*,2(v)\d ~ S*, 2 ((p)\d DR oc yjd - d DR , as seen in Fig. [31 

When considering dimensions smaller than dun, one 
must study the full dependence of the function 5k(p = 



<£> 2 /2, y). We show in Fig.@]the evolution of this function 
for p = 0. Starting from a constant function, one can 
clearly see that a linear cusp in y appears in a finite RG 
time |tz,|, close to 0.7 for the case shown. 

The appearance of a cusp along the flow leads to a 
breakdown of the superrotational invariance and of the 
associated WT identities. This is illustrated for the d = 3 
fixed point in Fig [5j there, z*(p) ^ S r ^(p) = 8*(p,0), 
which implies a breaking of the WT identity in Eq. (|27p . 
The different asymptotic behaviors at large p are eas- 
ily deduced from Eqs. (|35|36[) . from which we show that 
z*{p) ~ p-v/tt- 4 ^ and <5*(p,0) - p -{2v-v)/{d-i+fi) _ In 
d — 3, 2rj — fj is very small, which implies that the func- 
tion 8*(p,0) decreases slowly to 0. 

We also display in Figs. [6] and [7] the fixed-point solu- 
tions u'^(p) and z*(p) for different dimensions. One can 
see that the lower the dimension the steeper the curves. 
This means that for a numerical study of the critical 
properties in low dimensions, we need to discretize the 
field dependence in the NP-FRG equations with a small 
mesh, which entails a large number of points. The numer- 
ical integration is therefore more difficult and even be- 
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FIG. 4: NP-FRG flow of 5 k (0, y) for d = 4 < d Di j. The initial 
conditions at fc = A (i.e., t = 0) for u' k (p) and -Zfc(p) are taken 
at the fixed-point solution, u' t (p) and 2*(p), and that for dk 
is chosen as 5k=A(p,y) = 1- One observes that a linear cusp 
appears at a finite RG time. By construction 5k (0, 0) = I 
along the flow. 




FIG. 6: Fixed-point solution w»(p) for dimensions ranging 
from 3 (steepest) to 4.8 (smoother curve). The variable p 
has been rescaled by a factor e~ d so that all curves can be 
represented on the same plot. 
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P 

FIG. 5: Fixed-point solution in d = 3 for z,(p) (solid line) 
and 5*{p, 0) (dashed line). The two functions differ for a large 
enough field (by construction, they coincide at p = 0). 
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FIG. 7: Fixed-point solution z»(/o) for dimensions ranging 
from 3 (steepest) to 4.8 (smoother curve). The variable p 
has been rescaled by a factor e~ d so that all curves can be 
represented on the same plot. 



comes intractable in practice. Finally, we show in Fig. [5] 
the full dependence of the fixed-point solution S*(p,y). 

We now turn to the results concerning the critical ex- 
ponents. We begin with the anomalous dimensions 77 
and fj which are determined at the (critical) fixed point. 
(We recall that 77 characterizes the spatial decay of the 
"connected" pair correlation function and fj that of the 
"disconnected" pair correlation function at criticality.) 
Their dependence on the spatial dimension d is shown in 
Fig. [5] Above don — 5.1, we find that rj and fj rigor- 
ously coincide, fj = 77, which is the signature of dimen- 
sional reduction. (Their value is moreover equal to that 
found in the same order of the derivative expansion for 
the pure Ising model in dimension d — 2.) Below don, 
the two anomalous dimensions bifurcate and fj r\. Note 
that the predicted values satisfy the required bounds, i.e., 
277 > fj > r\, rj > (4 - d)/2, fj>A-d. 

We give in Table Q] our estimates for 77 and fj in d = 3 
and d = 4, where a comparison is possible with numerical 



determinations from computer studies (Monte Carlo cal- 
culations and T = ground-state determinations). One 
can see that the agreement is good or very good. (Real- 
space RG studies in d — 3^— provide results in the same 
range, with an anomalous dimension rj between 0.51 to 
0.56 and an exponent 8 between 1.50 and 1.56.) As is 
well known from the study of simpler models such as the 
pure O(N) model, the way to further improve the accu- 
racy of the exponents would be to consider higher orders 
of the approximation scheme presented in section IV-A. 

We have not been able to perform our NP-FRG cal- 
culation down to the lower critical dimension d = 2, 
where the values rj — 1 and fj = 2 are exactly known. 
The numerical resolution of the flow equations become 
extremely arduous in low dimension where the anoma- 
lous dimensions become large and approach their lower 
bound, (4 — d)/2 and 4 — d respectively. In addition, one 
also encounters numerical difficulties as one approaches 
the critical dimension dun from below. Indeed, as men- 
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FIG. 8: Fixed-point solution S,(p,y) in d = 3. Note that the 
point (p = 0, y = 0) corresponds to the lower right corner. 
The function is even in y (with a cusp around y = 0) and the 
part corresponding to y < is not shown. 



always valid. 

Finally, we also display our results in d = 3 and d = 4 
for the correlation length exponent v 1 which is associated 
with the relevant direction around the critical fixed point, 
in Table I. The agreement with the available data from 
computer studies is excellent in d = 4 and fair in d = 3. 
(Note that in the real-space RG studies in d = 3^—, the 
exponent v is in general too large, with values ranging 
from 1.39 to 2.25.) All the other critical exponents are 
obtained in the NP-FRG through the expected relations: 
7 = (2 - rj)u, 7 = (4 - fj)v, S={d + 4- fj)/(d - 4 + fj), 
etc^ 

The validity of the theoretical description is clearly 
confirmed by both the overall consistency and the quan- 
titative accuracy of the predictions in d — 3 and d = 4. 

B. Role of the auxiliary temperature and 
correction to "Grassmannian ultralocality" 



T), f] 




0.2 - 



FIG. 9: Anomalous dimensions rj and fj as functions of the 
spatial dimension d. The dashed lines are lower bounds for 
the anomalous dimensions [(4 — d)/2 and 4 — d for rj and fj 
respectively] and the crosses correspond to predictions from 
computer simulations and ground-state determinations. 



tioned in section IV-D, it is anticipated that the limit 
d — > dp R , y — > is nonuniform with a boundary layer 
in y 2 1 (don — d). More work will be needed to solve this 
boundary-layer problem numerically. 

Fig. [§] also provides evidence that the claim accord- 
ing to which the two exponents r\ and fj are related by a 
fixed ratio, fj = 2r/) ^ 24 ' 66 cannot be right in general. It 
is true that the relation fj = 2rj is exact in d — 2 (and 
according to a phenomenological RG, 67 also at first order 
in d = 2 + e) and is very closely obeyed by the numerical 
estimates in d — 3 and d — 4 (see Table H]). However, the 
latter type of "numerical evidence" can always be chal- 
lenged, and it is actually impossible to reach a definite 
conclusion on the sole basis of numerical results in se- 
lected dimensions (d — 3,4,5). On the other hand, the 
overall and continuous dependence on spatial dimension 
that we provide through the NP-FRG leads to a different 
and firmer type of answer: since one goes from fj — r/ for 
d > dun to fj > rj above, the relation fj — 2rj cannot be 



To assess the effect of a nonzero auxiliary tempera- 
ture and of a deviation from "Grassmannian ultralocal- 
ity", we have considered the simplest "non-ultralocal" 
contribution to the cumulants introduced in section VII- 
B of paper III. Combined with the above truncation of 
the "ultralocal" contributions of effective average action 
[Eqs. (|28l29p ]. this amounts to expressing the first cumu- 
lant as (see Eq. (136) of paper III) 



r fcl [#i] = jf (r*i [*!&)] + i(i + /?Mi)x 



Y k ($! (x, 00)0* £1)3* *i h) 



(52) 



with rfci[</>] given in Eq. (|28p. whereas the second-order 
cumulant is taken as purely "ultralocal" and given in 



d 


77 (this work) 


r\ (ground state) 


rj (Monte Carlo) 


3 


0.57(2) 


0.51 a 


0.5-0.6 / ' 9 ' h 


4 


0.24(1) 


0.18 b ,0.20 c ? 




d 


fj (this work) 


fj (ground state) 


fj (Monte Carlo) 


3 


1.08(2) 


1.02 a ,1.06 d ,l.l e 


1.0-1.04 / ' 9 '' 1 


4 


0.40(?) 


0.37 b ,0.38 c 




d 


v (this work) 


v (ground state) 


v (Monte Carlo) 


3 


1.08(2) 


1.0 e ,l.l\1.19 d 


1.0 m , 1.1 s , 1.31" 






,1.22 s ', 1.25 fc ,1.36 ! 








1.37 a 




4 


0.81(3) 


0.78 6 ,0.82 c 





TABLE I: Comparison between the anomalous dimensions 
obtained in the present work and in computer studies (a: 
b: El; c: H; d: JM; e: El; f: pi; g: [H; h: [13; i: 
j: [561 ]: k: [571 ]: 1: [581 ]: m: [5Sj]; n: [6(|). The error bars on the 
results from simulations and ground-state computations can 
be found in the original papers. The values of the temperature 
exponent 9 = 2 + rj — fj that we obtain are 1.49 in d = 3 and 
1.84 in d = 4, in excellent agreement with the other studies. 
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Eq. <J29J). The FRG flow equation for the function Y k (<f>) 
is directly obtained from Eq. (139) of paper III. In a 
dimensionless form, it reads 

d — 4 + f) k 

dtVk{<p) = ( 2 Vk - f)k)yk{<f>) h 7, t py'k( i p)+ 



fab) - ) {(d - 2)l<f- 2 \v) + 2[S k , (<p) 

- z k {v)]l<f{y) + (fj k - + 0(T k ), 

(53) 

where the threshold functions are defined in Appendix |B| 
In deriving the above expression from Eq. (139) of the 
preceding paper, we have used the dimensionless quantity 
y k defined through Y k (<f>) = (Z k k 2 /T k )y k (<p) and the fact 
that as k —> 0, 



} k (q 2 ,0)-P k (q 2 



1 



Y k (4>) Vk(<p)' 



(54) 



The above equation is valid provided that P^ which 
asymptotically goes as k 2 ~ n is subdominant with respect 
to Yfe. As discuss just below, this is indeed true, although 
y k goes to zero at the fixed point. 

The flow of the "ultralocal" quantities in the limit 
/3 _1 = has already been solved and we can introduce 
the corresponding solution in Eq. (|5"3")l . We have solved 
the latter equation and found that y k (<p) goes to zero as 
one approaches the fixed point, in such a way that the 
dimensionful quantity Y k (ip) goes to a constant Y*: this 
is illustrated in Fig. 10 for d = 3. 



Yu 




FIG. 10: Flow of the (dimensionful) function Y k of the di- 
mensionless field with p = ip 2 /2. After a transient regime, 
the function tends to a constant. 

In addition, we can investigate the effect of the "non- 
ultralocal" contribution to the flow of the "ultralocal" 
quantities when the auxiliary temperature is dif- 
ferent from zero. This can be done by considering the 
ERGE's for the "ultralocal" components of the first and 
the second cumulants given in Eqs. (134) and (135) of 
paper III. From the above results (in particular, the fact 
that Yfc goes to a finite constant when k Q), we know 



that Q k (q 2 , 4>) is asymptotically equal to P k (q 2 , <fi) : see 
Eq. and discussion below. The flow of the dimension- 
loss first cumulant therefore follows Eq. (151) of paper III 
and the contribution of the "non-ultralocal" piece goes 
as T k times a well-behaved function(al), with T k ~ k e / (3 
and 9 = 2 + i]-f]>0. 

The flow of the second cumulant is also of the form of 
Eq. (151) in the preceding paper, but the dimensionless 
beta function due to the "non-ultralocal" piece is now po- 
tentially singular when a cusp appears in the limit of zero 
(auxiliary) temperature of the second cumulant. More 
specifically, one finds that 



d t S k (ip\,V2 
yft / { 



Pi °°(<Pl,<P2)- 



<Pi 



z k (ipi)q 2 + s(q 2 ) + u'Ktp-i 



-perm(12)|, 
(55) 



where f3g °°(<pi, (f^) is the zero (auxiliary) temperature 
beta functional obtained by taking derivatives with re- 
spect to ipi and f 2 of the right-hand side of Eq. (|36|) 
and we have omitted subdominant terms coming from 
the difference between Q k and P k (see above). 

We can now follow the derivation done for the effect 
of the (bath) temperature in the formalism without su- 
perfields (see paper 11^). After changing variables from 
(pi,tp2 to x = (ipi + (p 2 )/2 and y=(<p 1 — ip 2 )/2, one ob- 
tains the solution for 5 k (x,y) in the close vicinity of the 
fixed point, when both T k and y go to zero, in the form 
of a boundary layer 



h(x,y) = 



A(x)- B{x)\ l + C{xY 



0(T 2 ), 
(56) 



where A(x), B(x) , C(x) are well-behaved functions of 
x and v4(a;)|C(x)| is equal to the absolute value of the 
coefficient of the cusp in the zero-temperature fixed-point 
function 8* (x, y). 

The connection between this boundary-layer phe- 
nomenon and the presence of rare low-energy excitations 
above the ground-state known as "droplets" 69 i 70 is the 
same as that already discussed in our previous work^ (see 
also Refs. [4.5]) and the discussion is not repeated here. 



VI. CONCLUSION 

In this paper and the preceding onei^ we have extended 
our NP-FRG approach of disordered systems, which was 
initiated in the previous articles of this series The ob- 
jective was to discuss the property of dimensional reduc- 
tion and its breakdown in the RFIM from the Parisi- 
Sourlas 1 ^ perspective of an underlying supersymmetry 
and its breaking. We have reformulated the NP-FRG 
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in a superfield formalism and, to do so, we have pro- 
posed a solution for properly selecting the ground state 
among the many metastable states that are present at 
zero temperature. Through the introduction of an appro- 
priate regulator and a supersymmetry-compatible non- 
perturbative approximation, we have been able to follow 
the supersymmetry (more precisely, the superrotational 
invariance) and its spontaneous breaking along the RG 
flow. 

Despite the fact that the effective hamiltonian (or mi- 
croscopic action) in the presence of a random field has 
numerous minima in the region of interest near the crit- 
ical point, dimensional reduction need not be systemat- 
ically broken. By implementing a NP-FRG flow for the 
cumulants of the renormalized disorder, our work shows 
that there is a finite range of dimension below the up- 
per critical dimension, duR, < d < d uc = 6, for which 
this multiplicity of minima has no effect on the long- 
distance properties of the model. More precisely, the 
scaling behavior around the critical point conforms to 
the dimensional-reduction predictions. The associated 
fixed point is characterized by a nonanalytic dependence 
of the effective action in the dimensionless fields, at odds 
with a naive description based on perturbation theory, 
but the nonanalyticity is too weak to alter the spectrum 
of critical exponents. It is only by going below a crit- 
ical dimension dun, which we numerically find around 
5, that spontaneous breaking of the supersymmetry (su- 
perrotational invariance) takes place. The robustness of 
our theoretical description, which explains the pending 
puzzles concerning the critical behavior of the RFIM, is 
supported by the good agreement obtained between the 
NP-FRG predictions for the critical exponents and the 
available values from computer studies. 

Finally, one may wonder whether the formalism devel- 
oped here can be useful in a different context. An ob- 
vious extension is a study of the long-distance physics 
of self-avoiding branched polymers and the associated 
property of reduction to the Yang-Lee edge singularity 
problem in two fewer dimensions. 71 Dimensional reduc- 
tion is known to hold in this casej 7 - 2 - and this provides 
a benchmark model to test the ability of our nonpertur- 
bative approximation scheme to reproduce this property. 
Another extension concerns the hysteresis behavior and 
out-of-equilibrium phase transitions of the RFIM at zero 
temperature when driven by a slow change of the exter- 
nal magnetic fieldi 7 -^— Work in those directions is under 
way. 

More challenging, and far more speculative too, is 
the question known as the "Gribov ambiguity" in the 
nonperturbative quantization of nonabelian gauge field 
theories! 26 ! 77 ' 78 The standard Faddeev-Popov gauge- 
fixing procedure aims at restricting the functional inte- 
gral over the gauge field to nonequivalent gauge config- 
urations (i.e, configurations that are not obtainable one 
from another by a gauge transformation). Unfortunately, 
the gauge-fixing conditions that respect Lorentz invari- 
ance and internal (color) symmetry, such as the Landau 



gauge condition, are inconsistent because they admit so- 
lutions that are equivalent up to a gauge transformation 
("Gribov copies"). There has been a strong theoretical 
effort to overcome this problem. This "Gribov ambigu- 
ity" does not affect the perturbative regime, so that cal- 
culations at high energy can be performed in the stan- 
dard Faddeev-Popov approach. However, very little is 
known on its influence in the nonperturbative regime. 
The Gribov-Zwanziger mode l 78 ' 79 enables one to reduce 
the number of copies taken into account in the functional, 
but not to a single one, and, up to now, there is no un- 
ambiguous gauge-fixing procedure. The strong connec- 
tions between the Faddeev-Popov and the Parisi-Sourlas 
formalisms, and their common failure when the relevant 
field equation has a multiplicity of solutions, have been 
underlined^ 6 - The relevance of the tools and concepts de- 
veloped in the present and the preceding papers to the 
problem of the Gribov ambiguity has yet to be investi- 
gated. 

Appendix A: Stability of the nonanalytic behavior 
along the RG flow 

To illustrate that a linear cusp in rj^. X2 (0 1 ' does 
not generate a supercusp and may break dimensional re- 
duction, we consider the ERGE for the Fourier transform 
of the 2-point proper vertex T^^iQ 2 ] 4>i,4>2) (for uniform 
fields) in the limit 02 — > <t>i- We switch to the variables 
4> = (0i + </>2)/2 and Y = (0i — <j>2)/2 and we study a 
situation in which a cusp is present: 

rf 2 V; 6 Y) = r^tf; ^ + rg> { q 2 -, <j>)\Y\ a + ■ ■ ■ , 

(A1) 

where Y — > and a < 2 (recall that is even in <fi 

and Y separately). With the assumption concerning the 
behavior of the proper vertices detailed in Sec. III-B, we 
find that the flow of , where for simplicity we take 
q 2 = 0, can be written in the Y — s- limit: 



dtT^ {q 2 = 0; (j), Y) ~ regular - 




J q' J q' 

- Ma"Wo(* - Y, 4> - y, <M\4*=4+y)] , 

(A2) 

where regular denotes the terms that could be obtained 
by directly considering the flow of rl^ (q 2 = 0; 4>, Y = 0) 
and by assuming a regular behavior of its beta-functional 
when Y — > 0. The contribution due to the derivatives 
of the third cumulant 1^3 is an even function of Y; 
the Y = term contributes to the regular component of 
the beta- functional and the next term is a 0(|F| Q ) with 
a' > 2. 

The above equation indicates that a "supercusp" with 
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a < 1 leads to an ill-defined flow for {q 2 — 0; 4>, Y = 
0). For a > 1 (cusp) one should also study the ERGE 
for the third cumulant to determine the leading nonana- 
lyticity (i.e., a') and therefore study the stability of the 
cusp under further RG flow. In any case, Eq. (|A2[) shows 
that a cusp weaker than linear, with 1 < a < 2, has no 
effect on the flow of rli 1 (q 2 = 0',(f>,Y = 0) and, since 
2??fc — rjk is obtained from the latter, is not expected to 
modify dimensional reduction. A linear cusp (a = 1) 
on the other hand provides a possible mechanism for the 
failure of dimensional reduction while not a priori gen- 
erating stronger nonanalyticities. In order to actually 
lead to breakdown of dimensional reduction, the linear 
cusp must of course remain in the renormalized second 
cumulant up to the appropriate fixed point. 



Appendix B: Dimensionless threshold functions 

The various dimensionless threshold functions corre- 
spond to the various 1-loop integrals involving the in- 
frared cutoff functions and the propagators (after account 
of the scaling dimensions) . More specifically, they are de- 
fined as followsi2&29^ 



dyy d l 2 - l d t 



91.92 
1 

2./o 



(pi(y) + w 1 )^{p 2 {y) + w 2 ) 92 



(Bl) 



compared to R when dimensional reduction is broken. 
These threshold functions then always appear multiplied 
by (77 — fj) in the flow equations. We define: 



1 dv .4/2-1 (r(y)+yr'(y)r 

10 {pi(y) + w 1 )^{ P2 {y) + w 2 )^ 

and 

h q d \w;z) = 

' , d/2 r (v) + y r '(y) f 2z + r (y) + v r '(y) 

„ dyV ( P {y) + W y I 2 
-^y— {z + r(y)+yr'(y))}. 



(B4) 



(B5) 



We summarize here a number of relations that 
are useful in the developments of Sec. IV. One 

has lq (w; z) — lqf] q2 (w,w; z, z), m q (w; z) = 
m qi ,q 2 (w,w;z,z), n q '(w;z) = n gi ,g 2 {w,w;z,z), 
j q a ' d \w;z) = j q ^ q l(w,w;z,z), with q 1 + q 2 = q. 
In addition, the functions rig 's and l q 's are related by 



n q d \w-,z) = ^- Y) l q a \{w-z). 



(d) 



(B6) 



From the above definitions, one also straightforwardly 
finds that 

M?!«s(i,2) = -?i&i, 92 (m). ( B7a ) 



m q d ] q2 (w 1 ,w 2 ;z 1 ,z 2 ) 

1 dyy^-% 



y{d y Pi(y)Y 



(pi(y) + w 1 )^{ P2 (y) + w 2 ) 92 j ' 

(B2) 



and 



d z J q fn,2) = - qi l q d Xl](l,2), (B7b) 



d Wl m q d ] q2 (l,2) - - qi m q d) +ltq2 (l,2), (B8a) 



(wi,w 2 ; zi,z 2 ) 



dyy d ' 2 - l d t 



ydyPiiy) 



(B3) 



with p a (y) = y(z a + r(y)),a = 1,2, y = q 2 /k 2 , and d t 
acts only on the dimensionless function r(y) — s(y)/y = 
Rk{yk 2 )/Zkyk 2 that is contained in the p a (y) ' a 27 i 28 i 40 (by 
definition, d t r{y) = —\qr{y) + 2yr'{y)\). Note that, in 
addition to the dependence onioE u'k(4>) and z = Zk{<f), 
the threshold functions explicitly depend on the scale k 
via the running anomalous dimension rjk- (Note also that 
the functions m^ d > and rS d ^ are not symmetric in the 
exchange of the indices 1 and 2.) The generic properties 
of these dimensionless threshold functions are discussed 
in detail in Ref. (HHSl. 

In the present study, it is also necessary to introduce 
additional threshold functions, which accounts for the 
fact that the cutoff function R has an anomalous scaling 



d Zl m q % 2 (l,2) = -qim q d Xl q2 (1,2) +—1^ (1,2), 

(B8b) 

d Z2 m q d ] q2 (l,2) = -^+^(1,2). (B8c) 

Similar relations hold for the functions Jg"'^, and h q d \ 
From these expressions one can obtain the derivatives 
with respect to the field arguments tp a by using that w = 
u'l(4>) and z = Zk{4>)- 



Appendix C: Some technical aspects of the 
numerical resolution 



1. Optimization of the IR cutoff function 

In the exact formulation of the NP-FRG, the results 
of the flow equations are independent of the IR regula- 
tor, which is only an intermediate means that does not 
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FIG. 11: Dependence of the anomalous dimensions rj (up- 
per panel), f\ (middle panel) and of the critical exponent v 
(lower panel) on the parameter a of the IR cutoff function in 
Eq. (|51[) . The "principle of minimal sensitivity" leads to a de- 
termination of the anomalous exponents r\ — 0.565, fj = 1.075, 
v = 1.08 (at the minima). The variations of the anomalous 
dimensions are very small over a wide range of a and enables 
us to evaluate the precision to A77 — 0.03 and Afj — 0.04. v 
is more sensitive to the regulator parameters and we estimate 
Av = §.2 



affect the long-distance physics. Once approximations 



are introduced, a residual dependence on the choice of 
regulator however remains. (A similar situation occurs 
in perturbation theory where a residual dependence on 
the parameters of the Borel resummation procedure is 
observed.) One can in some sense optimize the choice of 
IR cutoff function by demanding that the output, say the 
critical exponents, satisfies a property of "minimal sen- 
sitivity" such that by varying the characteristics of the 
cutoff function around the optimum, a minimal variation 
of the exponents results. This guarantees the stability of 
the results. We apply this procedure to the type of func- 
tion given in Eq. (|5ip where the parameters a, b, and c 
can be varied. 

We illustrate in Fig. [11] the dependence of the anoma- 
lous dimensions in d = 3 on the parameter a when b and 
c are kept fixed at the values 0.81 and 0.14 respectively 
(these latter values were determined by preliminary vari- 
ational studies to find minimal sensitivity). One can see 
that there is a large domain of values in which 77 and fj 
vary little (say by less than 5%). We use the "principle 
of minimal sensitivity" to determine the best estimate of 
the critical exponents [note that the minima of 77(a) and 
77(a) happen at very close values of a]. We also estimate 
the error bars from the observation that when a varies in 
the range [1.5,3], the maximum variation of 77 is of 0.03 
and that of 77 is of 0.04. These are the values that we 
have reported in Sec. [V] 



2. Regular izat ion of the numerical instabilities 

During the numerical integration of the flow, we en- 
countered numerical instabilities in physically unimpor- 
tant regions (typically, at large fields). The origin of these 
instabilities can be understood as follows. In the region 
of large fields, the threshold functions rapidly decrease 
to zero. In this limit, the flow of a function is given by 
the dimensional part, i.e. the flow obtained by setting 
v d ee [e.g. d t u' k (<p) ~ —\{d- 2r/ k + fj k )u' k (ip) + \{d - 
4+fjk)(f u'l((p)]. This flow no longer depends on the sec- 
ond derivative of the associated function with respect to 
the field. However, this diffusive-like term [u' k "(ip) in our 
previous example] is very important to stabilize the nu- 
merical integration. In practice, when necessary, we have 
therefore added by hand a small diffusive contribution to 
the flow equations and checked that the numerical re- 
sults are stable under varying the strength of these extra 
terms. 
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